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ABSTRACT 

Observational studies of nearby galaxies have demonstrated correlations between the mass of the 
central supermassive black holes (BHs) and properties of the host galaxies, notably the stellar bulge 
mass or central stellar velocity dispersion. Motivated by these correlations, the theoretical paradigm 
has emerged, in which BHs and bulges co-evolve. However, this picture was challenged by observational 
and theoretical studies, which hinted that the fundamental connection may be between BHs and dark 
matter halos, and not necessarily with their host galaxies. Based on a study of 3130 elliptical galaxies 
- selected from the Sloan Digital and ROSAT All Sky Surveys - we demonstrate that the central 
stellar velocity dispersion exhibits a significantly tighter correlation with the total gravitating mass, 
traced by the X-ray luminosity of the hot gas, than with the stellar mass. This hints that the central 
stellar velocity dispersion, and hence the central gravitational potential, may be the fundamental 
property of elliptical galaxies that is most tightly connected to the larger-scale dark matter halo. 
Furthermore, using the central stellar velocity dispersion as a surrogate for the BH mass, we find that 
in elliptical galaxies the inferred BH mass and inferred total gravitating mass within the virial radius 

1 g + 0 - 7 

(or within five effective radii) can be expressed as Mbh oc M tot -0,4 (or Mbh oc M 5rrf “ 0,6 ). These 
results are consistent with a picture in which the BH mass is directly set by the central stellar velocity 
dispersion, which, in turn, is determined by the total gravitating mass of the system. 

Subject headings: galaxies: elliptical and lenticular, cD, - galaxies: evolution, - galaxies: halos, - 
X-rays: galaxies, - X-rays: ISM 


1. INTRODUCTION 

Supermassive black holes (BHs), located at the center 
of every massive galaxy, are believed to have a profound 
effect on the evolution of their host galaxies. Rapidly 
growing BHs shine as active galactic nuclei (AGNs), and 
can release copious amounts of energy into their sur¬ 
roundings, which can heat and expel their gas supply, 
thereby quenching star formati o n and the BH growth 
dSilk fc Reesl ThMlKiiid 120031 iWvithe fc Loebl [20031 
Hopkin s et all I2006T) . Energetic feedback from BHs was 


suggested to play a pivotal role in buildin g the present- 
day observed galaxy lu minosity function dCroton et al.l 
120061 : iFaber et al .1120071 ). and hence BHs are a key com¬ 
ponent of modern galaxy formation theories. 

Observational studies of nearby galaxies have pointed 
out the existence of scaling relations between the BH 
mass (Mbh) and various properties of the host galax- 
ies, in particular with the stellar bulge ma ss, Mbulge 

(jMagorrian et, al.l 119081: lHaring fc Rlxl 20041) j_or cen- 

tral stellar velocity dispersion, a r dGebhar dt et, al.l[20001 
iFerrarese fc Merritt] 120001 : IGiiltekin et all 12009( 1 . To ex¬ 
plain the existence of these relatively tight scaling rela¬ 
tions, a widely accepted theoretical paradigm has been 
established, in which BHs grow in a self-regulated man¬ 
ner and co-evolve with the ir host bulges from early 
epochs to the present dav (ISilk fc Reesl 119981 iFabianl 
119991 iKingl 120031 iHoukins et all 120061) . However, sev¬ 
eral observational and theoretical studies revealed cor¬ 
relations between the BH mass and the total gravitat¬ 
ing mass, thereby hinti ng that dark matter halos may 
govern the BH growth (IFerra rese 2002; Bandara et all 


[20091 iBooth fc Schavel [2(TIol I20H [Bogdan et al.l 121)11 . 
Follow-up studies demonstrated that in late-type or 
bulgeless galaxies the dark m atter halo and t he BH 
mass do n ot corr elate tightly (IKormendv &; Benderll201ll 
ISun et aD 120131) . The absence of a tight relation in 
these galaxies may be due to low-mass or disk-dominated 
galaxies being less efficient at growing BHs, and dynami¬ 
cal effects that increase the scatter in th e Mr h — M to t re¬ 
lation for low-mass systems (|Volonteri et al.il2011[) . How¬ 
ever, observationally it is not well established, whether 
dark matter halos play a major role in driving the BH 
growth in elliptical galaxies. Therefore, to fully explore 
the connection between the BH mass and total halo mass, 
it is essential to perform a dedicated study focusing on 
elliptical galaxies. 

The main difficulty in conclusively determining 
whether the growth of BHs in elliptical galaxies is driven 
by stellar bulges or their large-scale dark matter halos 
arises from the demanding nature of BH and halo mass 
measurements. Specifically, dynamical BH mass mea¬ 
surements have only been obtained for ~ 80 galaxies, 
with abo ut half of these BHs being located in elliptica l 
galaxies (IGiiltekin et al.l 120091 : IMcConnell fe Mai I2013D . 
Of these syst ems, only a handful ha ve precisely measured 
halo masses (iHumphrev et al.]l2006D , which is insufficient 
to probe the importance of dark matter halos. Addition¬ 
ally, the local sample of galaxies with available BH and 
halo mass measurements does not offer a representative 
sample because galaxies residing in clusters are overrep¬ 
resented and include several brightest cluster galaxies. 
Therefore, instead of using individual galaxies, we may 
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opt to study the average properties of a large elliptical 
galaxy sample. 

In this paper, we probe, based on a statistically sig¬ 
nificant galaxy sample, whether BHs located in ellipti¬ 
cal galaxies exhibit a correlation with their total grav¬ 
itating mass. Because the direct determination of the 
BH mass (Mbh) and halo mass is not feasible for a 
large galaxy sample, we rely on proxies. For massive 
BHs, both observational and theoretical considerations 
suggest that the central stellar velocity dispersi on pro¬ 
vides a good measure of the BH mass (iGiiltcki n et al.1 
200l| IMcConnell fe Mai 120131) : for details see Section 
15.21 In the framework of the Sloan Digital Sky Survey 
(SDSS), a large number of nearby galaxies have been 
explored in a uniform manner, and SDSS observations 
provide the required accurate velocity dispersion mea¬ 
surements. For total halo mass measurements, the X-ray 
luminosity of the hot gas can be used as a ro bust tracer 
(iMathews et alJl2006t iKim fc Fabbianoll2013l ); for details 
see Section IQ To trace the hot gas component of ellip¬ 
ticals, we employ data from the ROSAT All Sky Survey 
(RASS). However, because of the relatively short individ¬ 
ual RASS observations, most passive elliptical galaxies 
beyond z ~ 0.01 remain undetected. Therefore, the X- 
ray properties of individual galaxies cannot be explored 
with RASS. However, when using the position informa¬ 
tion of a sufficiently large sample of SDSS galaxies, and 
stacking the individual X-ray observations, the resulting 
average galaxy population may show significant X-ray 
emission. Thus, to measure the X-ray luminosity, and 
hence the total halo mass of SDSS galaxies, we perform a 
stacking analysis. A stacking analysis has the additional 
advantage of probing the average properties of the galaxy 
population. Such a study naturally removes the distort¬ 
ing effects of outliers, as well as guards against biased 
selection techniques that are inherent in many studies 
focusing on small numbers of individual objects. 

The paper is structured as follows. In Section 2 we 
introduce and describe the sample of elliptical galaxies. 
In Section 3 we present the analysis of the ROSAT data, 
and in Section 4 we present our main result and show 
that the central stellar velocity dispersion of elliptical 
galaxies exhibits a very tight correlation with the X-ray 
luminosity of the hot gaseous component. In Section 5 we 
discuss that the central velocity dispersion and the gas 
X-ray luminosity can robustly trace the BH and total 
halo masses, and hence our results also imply that the 
BH mass tightly correlates with the total halo mass. We 
conclude in Section 6. 

2. SAMPLE SELECTION 

To construct a statistically significant galaxy sample, 
we rely on the SDSS Data Release 10 and the RASS 
data. First, we selected all SDSS galaxies that are within 
the redshift range of z = 0.01 — 0.05, have a stellar 
mass measurement, and have a stellar velocity dispersion 
measurement using the MPA-JHU analysis of the SDSS 
DRIP catalog (iKauffmann et aO2004lBrinchmann et al.l 
12004 ). The lower redshift boundary was introduced to 
exclude the most nearby galaxies, which given their prox¬ 
imity appear as the brightest on the sky (see Section 
ED. The higher redshift boundary is necessary to ensure 
the reliability of the central velocity dispersion measure¬ 
ments. Specifically, at z > 0.05 the Call triplet feature 


(8498A; 8542A; 8662A) is redshifted out of the observ¬ 
able range of SDSS spectroscopy (3900-9100A), thereby 
leading to larger uncertainties in the velocity dispersion 
measurement for galaxies beyond z > 0.05. 

In this work, we focus on the population of elliptical 
galaxies. A particular advantage of studying ellipticals 
is that for these systems the stellar mass is equivalent 
to the stellar bulge mass. Therefore, our results are not 
affected by uncertainties associated with bulge-disk de¬ 
composition, such as with studies of late-type galaxies. 
To study elliptical galaxies in the SDSS catalog, it is 
essential to robustly determine the morphology of the 
sample galaxies. 

Traditionally, elliptical galaxies were identified based 
on their locati on in the galaxy color-magnitude dia¬ 
gram dBell et al.ll2004) . However, recent results demon¬ 
strate that both early-type and late-type gala xies are dis¬ 
tribu ted across the entire color range fSchawinsk i et all 
l20ll . making it difficult to identify early-type galaxies 
solely using their galaxy colors. Therefore, the morphol¬ 
ogy of the sample galaxies must be determined by visual 
inspection on a case-by-case basis. To carry out this task 
for our large galaxy sample, we rel ied o n t he da ta pro¬ 
vided by the Galaxy Zoo project dLintott et al.l [2008). 
In the framework of Galaxy Zoo, citizen scientists visu¬ 
ally inspected the morphology of a vast number of SDSS 
galaxies, allowing us to unbiasedly select elliptical galax¬ 
ies in a galaxy-color-independent manner. To ensure that 
our sample only contains ellipticals, we included those 
systems that were deemed to be as elliptical galaxies by 
at least 70% of the citizen scientists. This definition re¬ 
sulted in a parent sample of 4592 elliptical galaxies that 
are observed with SDSS and are within the redshift range 
of z = 0.01 — 0.05. Additionally, we visually inspected 
all 4592 galaxies, and in agreement with the Galaxy Zoo 
project found that all selected galaxies exhibit an ellip¬ 
tical morphology. As a caveat, we note that the sample 
may contain SO galaxies in projection, which are diffi¬ 
cult to remove robustly. Moreover, we established that 
none of the galaxies in our parent sample appears to be 
substantially star forming. 

According to previous observational studies, low- 
stellar-mass systems cannot retain a significant amount 
of h ot X-ray gas in th e ir gravitational potentia l 
well (IQ’Sullivan et al.l 120011 IBoedan fc Gilfanovl I2011D . 
Therefore, we excluded all galaxies with stellar masses 
below logM*/M 0 = 10.0. 

To ensure that any observed X-ray emission within the 
stacked data is not being produced by AGNs, we used 
two methods to identify and remove AGNs. First, we 
applied typical optical emission line diagnostic diagrams 
from the SDSS optical spectroscopy. Here we choose to 
use both the [Nii1/Ha-[Ollll/H/3 and [Oil/Ho-[()[ll|/Hb 
diagnostics dKewlev et al.ll200ra iSchawinski et, al.ll2007l) . 
When combined, these ionization ratio diagrams effi¬ 
ciently and cleanly separate AGN emission from lower- 
ionization emission, such as that from star formation. 
Those sources that lie above the relevant demarcation 
curves are identified as AGNs and removed. Second, all 
galaxies that were directly detected in the RASS im¬ 
ages were further removed from the sample. However, 
as evidenced by inspection of the SDSS optical images, 
these particular sources were the brightest cluster galax- 
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TABLE 1 

Observed Galaxy Properties in Each Bin. 


Mi, 

(10 10 M 0 ) 

(1) 

Mi, 

(10 10 M 0 ) 
(2) 

Co 

(km s -1 ) 

(3) 

Galaxy 

# 

(4) 

Distance 

(Mpc) 

(5) 

N h 

(10 20 cm" 2 ) 
(6) 

Exposure 

(ks) 

(7) 

tx.LMXB 
(10 40 erg s _1 ) 
(8) 

T'Xjgas 

(10 40 erg s _1 ) 

(9) 

1.00 - 2.51 

1.58 

78.7 

265 

168.9 

2.03 

117.2 

0.04 

- 

1.00 - 2.51 

1.70 

111.8 

501 

165.4 

2.04 

212.2 

0.04 

- 

1.00 - 2.51 

1.62 

145.5 

192 

169.1 

2.05 

85.8 

0.04 

- 

2.51 - 6.31 

3.89 

115.2 

246 

179.5 

2.07 

105.2 

0.09 

< 1.0 

2.51 - 6.31 

3.98 

149.1 

717 

172.0 

2.06 

301.1 

0.09 

1.2 ±0.8 

2.51 - 6.31 

3.80 

193.5 

321 

164.7 

2.03 

142.8 

0.09 

4.2 ± 1.1 

6.31 - 15.85 

9.12 

155.6 

102 

176.7 

2.08 

40.9 

0.20 

2.6 ± 2.2 

6.31 - 15.85 

9.12 

197.3 

567 

172.6 

2.04 

239.7 

0.22 

6.4 ±0.9 

6.31 - 15.85 

8.70 

254.0 

117 

161.9 

2.04 

48.2 

0.21 

12 ±2 

15.85 - 39.81 

19.50 

249.6 

102 

169.8 

2.07 

43.7 

0.47 

14 ±2 


Note. Columns are as follows. (1) Stellar mass range of galaxies in units of 10 10 Mq. (2) Median stellar mass in units of 10 10 Mq. Note 
that for ellipticals the stellar mass is equivalent with the stellar bulge mass. (3) Median central stellar velocity dispersion given by SDSS. 
(4) Number of g alaxies in each bin. (5) Median distance. (6) Median line-of-sight column density derived from the HI map provided by 
the LAB survey dKalberla et al.l|2005l 1. (7) Combined exposure time of the galaxies in each bin. (8) X-ray luminosity of the population of 
unresolved LMXBs in units of 10 40 erg s —1 . Their luminosity is derived from the average LMXB luminosity function (IGilfanovll20041) . (9) 
Observed X-ray luminosity of the hot gas in units of 10 40 erg s —1 . The contribution of LMXBs from Column (8) is subtracted. Note that 
Lx,gas for the first three rows are consistent with Lx,gas — 0- 


ies at (or close to) the center of large groups and clusters. 
Based on these methods, we identified and removed 527 
AGNs from the sample. 

In addition to these, we searched through the SDSS 
and RASS images and filtered 299 sources that were con¬ 
taminated by foreground or background objects, such as 
a nearby galaxy in projection or a luminous background 
AGN. Such sources could contaminate the stacks and 
hence must be removed. We also searched for galaxies 
that had a companion with similar redshift within a 3' 
projected distance. The dark matter halos of these inter¬ 
acting galaxies may not be virialized, and hence the ob¬ 
served X-ray luminosity cannot be robustly used to trace 
their halo mass. We identified 185 interacting sources in 
our sample. Finally, we removed those 451 galaxies that 
were located in rich groups or clusters. The detected X- 
ray emission in these galaxies would likely be dominated 
by the large-scale hot diffuse group/cluster gas and not 
by hot gaseous coronae of the sample galaxies. The final 
sample contains 3130 elliptical galaxies. 

3. ANALYSIS OF THE ROSAT DATA 

The RASS explored the entire sky in the soft X-ray 
regime, providing the only all-sky survey at X-ray ener¬ 
gies to date. In the framework of the RASS, 1378 fields 
were observed, the majority performed in scanning mode, 
which resulted in fields with 6.4° x 6.4° length sides. The 
average exposure time of the RASS observations is ap¬ 
proximately 430 s, and hence most non-AGN galaxies 
beyond z ~ 0.01 are not detected. However, despite 
the shallow average exposure times of the RASS observa¬ 
tions, the most nearby galaxies may be detected. These 
individual detections could bias the stacked signal, and 
hence they were excluded from our study. Because the 
individual observations may have significantly shorter or 
longer exposures than the average value, it is necessary 
to account for the variations in the exposure times. To 
this end, we utilized the exposure maps that are provided 
with each observation. The combined exposure time of 
our 3130 sample galaxies amounts to ~ 1.3 Msec. 


Although ROSAT collected photons across the 0.1—2.4 
keV energy range, in our analysis we opted to use the 
0.4 — 2.4 keV band for the following three reasons. First, 
the ROSAT effective area sharply decreases below 0.4 
keV, resulting in fewer detected source photons at lower 
energies. Second, the instrumental background level is 
relatively high in the 0.1 — 0.4 keV band, so the signal-to- 
noise ratio may be reduced when the softest energy range 
is also used. Third, given the average level of Galactic 
absorption (IVh = 2.1 x 10 20 cm~ 2 ) for our galaxy sam¬ 
ple, we do not expect significant emission from the hot 
gaseous component below 0.4 keV energy. 

Compared to present-day X-ray telescopes, the point- 
spread function of ROSAT is relatively broad, which 
must be considered when choosing the source and back¬ 
ground extraction regions. To this end, we used circular 
regions with a four-pixel (180") radius centered on the 
stacked images to extract the source region. This region 
confines about 83% o f the source counts dBoesel [200(1 
iBohringer et al.ll2014D . This region is optimal because 
it contains a large fraction of the source counts, while 
it also minimizes the number of background counts. We 
note that choosing a somewhat smaller or larger source 
extraction region does not affect our results in any sig¬ 
nificant way. To account for the source counts falling 
outside of our extraction region, we increased the ob¬ 
served net count rates by 1/0.83 = 1.20. We subtract 
the background components by using local background 
regions. Utilizing local background regions is advanta¬ 
geous because both the instrumental and sky background 
components can be precisely subtracted. Therefore, we 
used circular annuli with a 8 — 14 pixel (360" — 630") 
radii to account for the background components. Within 
these regions we do not expect a notable contribution 
from the stacked signal, thereby ensuring a precise back¬ 
ground subtraction. 

To perform our stacking analysis, we produced 24' x 24' 
cutout images of the individual X-ray images and the 
corresponding exposure maps. The cutout X-ray images 
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Fig. 1.— Analyzed sample of elliptical galaxies. Black dots show 
the velocity dispersion as a function of the stellar mass for the par¬ 
ent sample of 4592 elliptical galaxies observed by SDSS. The red 
dots depict the 3130 galaxies selected in our stacking analysis after 
removal of all interacting systems, galaxies located in rich groups 
or clusters, images contaminated by foreground or background ob¬ 
jects, and AGN detections. The dashed regions show the selection 
regions as defined by the stellar mass and central velocity disper¬ 
sion. 

(and their respective exposure maps) were summed to 
produce stacked X-ray images of the galaxies in our sam¬ 
ple for a given stellar mass and central stellar velocity 
dispersion (see Section 1X11) . 

4. RESULTS 

4.1. Binning the Galaxy Sample 

Our large galaxy sample allows us to probe whether 
the central stellar velocity dispersion and the X-ray lu¬ 
minosity exhibit a tight correlation with each other. To 
this end, we grouped the galaxy sample based on their 
stellar mass and central velocity dispersion. We split 
the sample into four equal-width stellar mass bins in 
the logM*/M 0 = 10.0 — 11.6 stellar mass range. The 
galaxies within each of these bins are further subdivided 
according to their velocity dispersion (Figure Q]). Specif¬ 
ically, we grouped 68% of galaxies that follow the av¬ 
erage stellar velocity dispersion-stellar mass relation in 
the central region of each mass bin, while galaxies with 
seemingly low or high velocity dispersions, relative to 
their stellar mass, are located below or above the cen¬ 
tral regions, respectively. Because of the low number of 
galaxies in the logM*/M 0 = 11.2 — 11.6 stellar mass 
range, we only defined one group. The average proper¬ 
ties of the galaxies residing in each bin are tabulated in 
Table □ 

Given that we use the stellar mass as a major crite¬ 
rion to divide the galaxy sample into multiple bins, we 
investigate the stellar mass distributions within each bin. 
The obtained stellar mass distributions are shown in Fig¬ 
ure [2l which demonstrates that within the same stellar 
mass interval, the stellar mass distributions are nearly 


10.0 10.1 10.2 10.3 10.4 10.4 10.5 10.6 10.7 10.8 



10.8 10.9 11.0 11.1 11.2 11.2 11.3 11.4 11.5 11.6 

Stellar Mass (log M/M s ) 


Fig. 2.— Four panels show the stellar mass distributions for each 
galaxy bin. Within each stellar mass bin the mass distributions 
are similar, independent of the median velocity dispersions. Based 
on Kolmogorov-Smirnov tests, we find no statistically significant 
evidence that any of the stellar velocity dispersion samples are 
drawn from different distributions. Note that in the log M*/Mg = 
11.2 — 11.6 stellar mass range, only one group is defined because of 
the low number of galaxies. 

identical. To statistically probe if the galaxy samples 
are drawn from the same parent population within each 
stellar mass bin, we performed a series of Kolmogorov- 
Smirnov tests. We obtained probabilities > 0.09 for the 
comparison of each low-, central-, and high-velocity dis¬ 
persion samples within the same stellar mass bins. This 
implies that we have no statistically significant evidence 
to suggest that any of the velocity dispersion samples are 
drawn from separate distributions. 

The average redshift of the galaxies in our sample is 
z ~ 0.0415, which corresponds to a distance of 170.1 
Mpc. To ensure that possible differences in the red- 
shift distribution are not affecting our results, we in¬ 
vestigated the distance distributions of each galaxy bin. 
In particular, we performed Kolmogorov-Smirnov tests, 
which pointed out that there is no evidence at the 
> 1% level that the log M*/M 0 = 10.0 — 10.4 and 
logM*/M 0 = 10.4 - 10.8 and logM*/M 0 = 10.4 - 10.8 
and logM*/M 0 = 10.8 — 11.2 mass bins are not drawn 
from the same distributions. However, statistically we 
would reject Hq at the 1% level for the comparison be¬ 
tween the log M*/M 0 = 10.0 — 10.4 and logM*/M 0 = 
10.8 — 11.2 samples because of a slightly larger popula¬ 
tion of more massive systems at higher redshift, which 
is to be expected for any flux-limited sample. But we 
would still accept Hq at the 0.1% level for these bins, 
so any difference can still be considered marginal. The 
median distances of the galaxies in each bin are close 
to 170.1 Mpc, and none of the bins deviates more than 
~ 5% from this average value (Table [T| . 

4.2. A Tight a c — Lx Correlation 
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Fig. 3.— Schematic diagram of stacked RASS X-ray images for the 10 bins defined by their stellar mass and velocity dispersions. No 
exposure correction is applied. The images were smoothed with a Gaussian kernel of three-pixel size. The circular region has a radius of 
four pixels and represents the source extraction region. The bins with high (> 190 km s" 1 ) median velocity dispersion exhibit a statistically 
significant, 4.2a — 7.3a, detection, and bins with low (< 190 km s _1 ) median velocity dispersion either show a moderate detection (~ 2a) 
or a nondetection. Note that in the lowest stellar mass bin (logM*/Mg = 10.0 — 10.4), no statistically significant source is detected, 
which is partly due to the shallow potential well of these galaxies, and partly due to the low effective depth of the stacked data in the bin 
characterized with the highest (149 km s —x ) velocity dispersion. 


For each galaxy bin, we stacked the X-ray photons in 
the 0.4 — 2.4 keV band RASS cutout images and their 
corresponding exposure maps. The stacked images are 
shown in Figure [3] which reveal statistically significant 
detections in several bins irrespective of the median stel¬ 
lar mass. In particular, within the stacked images statis¬ 
tically significant (4.2cr—7.3<r) X-ray detections are found 
for galaxies characterized with > 190 km s -1 median ve¬ 
locity dispersions. However, only moderate (~ 2<r) de¬ 
tections or nondetections are observed for galaxies ex¬ 
hibiting < 190 km s -1 median velocity dispersions. 

To convert the net count rates to flux, we assumed 
an optically thin thermal plasma emission model with 
kT = 0.7 keV and 0.3 solar abundance. Although this 
gas temperature and metal ab un dance are typical for 
elliptical galaxies (|Humohrev fe Buotell2006D . we inves¬ 
tigate whether deviations from these parameters could 
have any significant effect on our results. The gas 
temperature for galaxies in the stellar mass range of 
log A fi /Mg, = 10. 4— 11.6 is in the range of kT ~ 0.4 —0.9 
keV flBoedan fe Gilfanovll201ll) . Assuming that the lower 


or upper bounds of this temperature range are at fixed 
metal abundances would result in a ~ 3% difference 
in the counts-to-flux conversion. Whereas most ellip¬ 
ticals exhibit subsolar metal abundances, some massive 
galaxies (logM*/MQ > 11.0) in our sample may have 
abundances close to the solar value. Using solar abun¬ 
dances instead of a 0.3 solar value at a fixed temperature 
would result in a ~ 7% difference in the counts-to-flux 
conversion. To probe the combined effects of varying 
temperatures and abundances, we compared two mod¬ 
els: one with kT = 0.4 and 0.2 solar abundances, and 
another with kT = 0.9 and solar abundance. We found 
that the counts-to-flux conversions for these two extreme 
cases are only ~ 10% different. Therefore, our assumed 
model with kT = 0.7 keV and 0.3 solar abundance en¬ 
sures that our flux conversion is accurate to within ~ 6%. 
This value is lower than the statistical uncertainties, and 
hence does not influence our results in any significant 
way. Additionally, we note that the X-ray hardness ra¬ 
tios also suggest that the stacked signal is dominated 
by hot gas with temperatures of kT ~ 0.7 — 0.8 keV 
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Fig. 4.— Results of the stacking analysis. The left panel shows the median central stellar velocity dispersion (cr c ) as a function of the 
median stellar mass. The right panel depicts cr c as a function of the X-ray luminosity and reveals a near-perfect correlation. Note that 
the lowest stellar mass bin is not shown because no significant hot X-ray emitting gas is detected. The uncertainties in cr c represent the 
scatter (67th percentile) in the sample galaxies. The uncertainties in the X-ray luminosity are derived from 10 4 jackknife simulations of 
the galaxies within each bin. The depicted uncertainties are computed from the 67th percentile of the simulations. 


(Section 14.31) . The average X-ray luminosities were cal¬ 
culated assuming the median redshift of each bin. These 
luminosities were corrected for the median line-of-siglit 
column density for each bin (Table (TJ, which values were 
derived from the HI map provided by the LAB survey 
dKalberla et alJl2005lb 

To derive the X-ray luminosity associated with the hot 
gaseous emission, it is essential to subtract the contribu¬ 
tion of other X-ray emitting sources. The largest con¬ 
tribution comes from the population of low-mass X-ray 
binaries (LMXBs). In elliptical galaxies, the emission 
from LMX Bs is propor tional to the stellar mass of the 
galaxies (lGilfanovll2004f) . Therefore, we used the average 
LMXB luminosity function, the average LMXB spectrum 
(llrwin et al.l 120031 ). and the stellar mass of the galaxies 
in each bin to derive the X-ray flux originating from un¬ 
resolved LMXBs. We find that the population of unre¬ 
solved LMXBs plays a minor role in the stacked signal 
because in those bins where hot gas is detected their 
contribution remains at the 3 — 15% level. The emission 
from fainter compact objects, such as active binaries and 
cataclysmic variables, contributes a t the < 1 % level for 
the el liptical galaxies in our samp le (Bogdan fc Gilfanovl 
ML 12011 |s azonov et al.l I2006f) . Hence, their contri¬ 
bution to the stacked signal in our sample galaxies is 
virtually negligible. 

To derive the uncertainties in the X-ray luminosity, we 
utilized Monte Carlo (MC) random simulations. For each 
MC simulation, we randomly stacked the X-ray images 
of 90% of the sources present in each bin to construct the 
X-ray luminosity distribution of the stacked emission. 
From the 10 4 simulations, we assessed the median and 
the 67th percentiles of the distributions to measure the 
uncertainties in the stacks. This process guards against 


the contribution from extremely luminous and underlu- 
rninous sources in the sample that may dominate the 
observed X-ray flux. Additionally, we also derived the 
uncertainties using bootstrap analysis, which resulted in 
uncertainties identical to those obtained by the MC ran¬ 
dom simulations. 

The result of our stacking analysis is presented in Fig¬ 
ure [H which reveals that, for a given stellar mass, el¬ 
lipticals with the largest velocity dispersions, present the 
strongest X-ray luminosities. The best-fit relation be¬ 
tween the observed X-ray luminosity and central velocity 
dispersion can be described as: 

logLx(erg s _1 ) = (33.2±1.6) + (3.3±0.7) logcr c (km s _1 ) . 

Our X-ray stacking analysis demonstrates a near-perfect 
cr c — Lx correlation with a Pearson correlation coefficient 
of 0.98 and y 2 = 1.3 for five degrees of freedom. The 
critical value to reject Hq (i.e. no correlation between 
these quantities) at the 90% level is x 2 = 1.6, and hence 
we conclude at the ~ 90 — 95% confidence level that 
the central velocity dispersion and the X-ray luminosity 
tightly correlate with each other. 

4.3. X-ray Hardness Ratios 

To confirm that the detected X-ray signal originates 
from hot gaseous emission with sub-keV temperatures, 
we derive X-ray hardness ratios using the stacked RASS 
data. We define the hardness ratio as the ratio of count 
rates observed in the 1.0 — 2.4 keV and 0.4 — 1.0 keV 
energy ranges, that is HR = Ci.o- 2 . 4 kev/C 0 . 4 -i.okeV- To 
obtain statistically meaningful hardness ratios, we mea¬ 
sure them only in the largest median velocity dispersion 
(> 190 km s _1 ) bins. The hardness ratios are in the 
range of HR = 0.33 — 0.55, where the largest value corre- 
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sponds to more massive halos. These values are consis¬ 
tent with an emission spectrum dominated by hot ion¬ 
ized gas with kT ~ 0.6 — 0.8 keV temperature, for which 
HR = 0.35 — 0.5 is expected. Moreover, observations of 
nearby galaxies demonstrate that these gas temperatures 
are typica l for galax ies w ith stellar mass of ~ 10 11 M 0 
dBogidkn fc Gilfanovll2011lh 

The derived hardness ratios are in conflict with a model 
in which the observed X-ray emission originates from 
AGNs or unresolved X-ray binaries. Indeed, a power- 
law model with a slope of T = 1.7 , which is typical 
for AGNs or X-ray binaries, would result in HR « 1 
at these energy ranges. We also note that an obscured 
AGN would have HR > 1, which is also inconsistent with 
our observations (see Section l4~4l for details). Thus, the 
hardness ratios demonstrate that AGNs or other hard X- 
ray emitting components do not play a major role in the 
observed signal, but the observed X-ray emission is pre¬ 
sumably dominated by hot X-ray gas with kT ~ 0.6 —0.8 
keV temperature. 

4.4. Obscured AGNs are Not the Source of the 
Stacked X-ray Signal 

Because X-ray stacking has been routinely used to in¬ 
fer the presence of heavily obscured AGNs buried at 
the centers of distant galaxy populations, we investi¬ 
gate whether the X-ray signals identified in the stacked 
RASS images are arising due to the presence of ob¬ 
scured AGNs. Therefore, we make the opposing ansatz 
that the observed X-ray luminosity is due to obscured 
BH growth. Within our stacked images, we resolve 
an X-ray luminosity of ~ 3 x 10 4 ° erg s _1 . To con¬ 
vert this value to a bolometric X-ray luminosity, we as- 
su me a conversio n factor kbo\ = Lboi/Ao. 5 - 2 keV = 20 
(|Lusso et al.H2012f h and a power-law model with a slope 
of T = 1.9 and intrinsic Nr = 10 22 cm -2 , which is typi¬ 
cal of a mildly obscured Type-2 AGN (e.g. IRisaliti et al.1 
IMl). Thus, an observed (absorbed) X-ray luminosity 
of 3 x 10 40 erg s _1 corresponds to the AGN bolomet¬ 
ric luminosity of 1.2 x 10 42 erg s~ 4 , for those galaxies 
with BH masses of ~ 3 x 10 s Mq. Such a BH mass cor¬ 
responds to an Eddington limited accretion luminosity 
of ~ 4 x 10 46 erg s _1 , and hence, based on the observed 
X-ray luminosity, a lower limit of ~ 1CP 4 for the Edding¬ 
ton ratio. Such an Eddington ratio is typical of AGNs 
present in star-forming spiral galaxies, which have a plen¬ 
tiful supply of cool gas from which tcyform a radiatively 
efficient accretion disk (|Goulding et al.ll2010f) . Even pow¬ 
erful radio AGNs that are observed in elliptical galaxies 
are typica lly found to have far lower Eddington ratios 
: I In 2002 ). For AGN with Eddington ratios > 10 -4 , we 
would expect to readily observe optical AGNs emission 
lines. However, sources such as these were removed dur¬ 
ing our galaxy selection process. 

Furthermore, our X-ray images are constructed in the 
soft X-ray regime, E < 2.4 keV; such X-ray photons 
are particularly susceptible to even low levels of gas ab¬ 
sorption. Hence, for an AGN to be missed at optical 
wavelengths, this would require relatively large levels of 
extinction toward the central AGN, and any soft X-ray 
photons arising from this AGN would still be obscured in 
an X-ray stack. Indeed, all previous X-ray stacking anal¬ 
yses have harnessed only hard E > 4 keV photons to infer 


the presence of an obscured AGN popul ation, as the sof t 
photons remain hidden in the stacks dDaddi et al.ll2007f) . 
Thus, we conclude that the X-ray emission observed in 
the stacks does not arise because of the presence of ob¬ 
scured AGNs, and must instead originate from the pres¬ 
ence of galaxy-wide hot gas. 

4.5. Robustness of the Background Subtraction 

The nondetections of galaxies in the lowest stellar mass 
bin (log M s tai/M 0 = 10.0 — 10.4) demonstrate that our 
background subtraction is robust. Indeed, most galaxies 
within this stellar mass bin are not expected to reside 
in a sufficiently massive halo that is capableyof retainin g 
significant quantities of hot X-ray gas (lO’Sull ivan et ahl 
120011 : [Mathews et al.ll2006t iKim fc Fabbianol I2013T) . 

To further confirm that the detected signal is not of 
spurious origin, we also stacked random fields within the 
SDSS field of view. We selected 500 random coordi¬ 
nates covered by the SDSS, and stacked the correspond¬ 
ing RASS images following our previous procedure. The 
obtained stacked image of random fields consists of a 
flat background and lacks any evidence for a statistically 
significant detection at the center. Hence, the detected 
sources, associated with the position of the stacked galax¬ 
ies have a physical origin. 

5. DISCUSSION 

5.1. Central Stellar Velocity Dispersion Correlates with 

the Total Halo Mass 

X-ray observations have demonstrated that ellipti¬ 
cal galaxies with very similar optical properties ca n 
have v astly different X - ray p r operties dForman et al.1 
19851 [O ’Sulliva n et all feOOH IMathews et all 120031 
Bogdan fe Gilfanovl 1201111 . Specifically, these studies 
showed that the luminosity of the hot X-ray gas, which 
dominates the X-ray appearance of massive ellipticals, 
can vary by up to two orders of magnitude at a fixed 
stellar mass. However, the amount of hot X-ray gas 
that can be retained in the galaxy’s gravitational po¬ 
tential well is primarily determined by the dark matter 
halo mass rather than by the stellar mass. Motivated 
by this realization, several studies explored the connec¬ 
tion bet ween the X-ray luminosity and halo mass of e l- 
lipticals (IMathews et all 120061 : IKim fe Fabbianol 1201311 . 
These works demonstrated that higher X-ray luminosi¬ 
ties correspond to more massive halos, and hence the 
X-ray luminosity of the hot gas correlates with the halo 
mass. 

Indeed, IMathews et all (120061) established the scaling 
relation between the X-ray luminosity and virial mass 
for objects within the log Mtot/M© = 12.8 — 14.5 range, 
which sample included galaxie s, galaxy groups , and p oor 
galaxy clusters. In their study. IMathews et all (120061) re¬ 
lied on X-ray observations of the hot gas to derive the to¬ 
tal halo mass by assuming hydrostatic equilibrium. The 
obtained best-fit relation between the X-ray luminosity 
and the virial mass is described a s Lx oc M t 2 0 4 . In a 
follow-up study, IKim fe Fabbianol (1201311 investigated a 
sample of nearby early-type galaxies with a partic ula r fo¬ 
cus on gas-rich galaxies. In t heir work. IKim fc Fabbianol 
(|2013f) utilized the results of iDeason et al.1 (l2012f) . who 
determined the total mass within 5r e ff using optical kine¬ 
matics data of globular clusters and planetary nebulae. 
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Fig. 5.— Scaling relations between the inferred BH mass and the inferred total mass within the virial radius (left panel) or within 5r e ff 
(right panel). The best-fit Mbh — Mtot and Mbh — M 5 r eff relations are shown with the dashed lines. The uncertainties in total and halo 
mass are calculated by summing in quadrature the unc ertainties measured fr om t he MC simulations of th e X-ray luminosities, shown in 
Figure[4] with the rms scatter observed in the relations of I Mathews et ahl (1200611 and I Kim &; Fabbiancl (120131) . The uncertainties in Mbh are 
ca lculated by summing i n quadrature the 67th percentile scatter in our defined velocity dispersion bins with the observed Mbh — 0 c relation 
of IGiiltekin et ahl (|20091 1. The insets show the probability distributions of the slopes obtained from the fitting procedure. Overplotted are 
individual galaxies with dynamically measured BH mass and precisely measured total halo masses. 


iKim fe Fabbiand (|2013T> found that for gas-rich early- 
type galaxies with Lx, gas > 10 40 erg s _1 , which are rele¬ 
vant for our study, the Lx — -^ 5 r eff relation is extremely 
tight and its rms deviation is only 0.128 dex. Their best- 
fit relation can be expressed as Lx oc Mf r 3 . 

The existence of a tight Lx — M to t relation for grav¬ 
itationally collapsed struct ures is expect ed from simple 
self-similar considerations (iKaised 11986 *). Because the 
self-similarity applies not only for the collision less dark 
matter particles but also for the weakly collisional hot 
ioniz ed gas that is heated through shocks dNavarro et al.l 
we expect to observe power-law scaling relations 
between the X-ray properties of the gas. Specifically, 
under the assumption of hydrostatic equilibrium, the X- 
ray luminosity is predicted to trace the virial mass as 

Lx oc M 4 / 3 . While the observed relations are steeper 
than that expected from analytical models, this simple 
model illustrates that the observed correlations are phys¬ 
ically motivated. 

Based on these considerations we conclude that, on av¬ 
erage, the gas X-ray luminosity of our elliptical galaxy 
sample can be robustly used to trace the total halo mass. 
Combining this fact with the observed tight a c — Lx re¬ 
lation (Section 14.211 . we infer that for massive elliptical 
galaxies the central stellar velocity dispersion tightly cor¬ 
relates with the total halo mass. By definition, a c pro¬ 
vides a dynamical estimate of the central gravitational 
potential of the galaxy. As such, our results suggest that 
the depth of the central potential well and the total grav¬ 
itating mass of elliptical galaxies are intimately linked. 
Moreover, for elliptical galaxies, the er c — Lx relation 
is significantly tighter than the Mb u i ge ~ Ax correlation 
(Section [53]). This implies that the central stellar veloc¬ 
ity dispersion may be the fundamental galaxy property 
that is most tightly connected to the total halo mass. 


This in turn hints that for elliptical galaxies a c can more 
accurately trace the dark matter halos than the stellar 
mass. 

5.2. BH Mass Correlates with the Total Halo Mass 

Observational studies have repeatedly pointed out 
that the central stellar velocity dispersion correlates 
closely with the BH mass, and hence a r is a widely 
used surrogate for the BH mass (Gebh ardt et al.l 120001 : 

I Tremaine et al.l l2002t iKormendv fe Bended 1201 ill . Re¬ 
cent studies employing larger galaxy samples demon¬ 
strated that the Mbh — cr c scalin g relati on is particu¬ 
larly tight for elliptica l galaxies (IGiiltekin et al.1 120091 : 
iMcGonnell fc Mall2013D , which are relevant for our study. 
In these works the central velocity dispersion refers to 
the luminosity-weighted line-of-sight velocity dispersion 
measured within the half-light radius of the galaxy. The 
velocity dispersions provided by SDSS are measured 
within the 3" spectroscopic fiber. Our sample galaxies lie 
at a median distance of D = 170 Mpc, where the SDSS 
fiber size corresponds to a projected radius of ~ 1.2 kpc. 
The typical half-light radius of elliptical galaxies in the 
studied stellar mass range is a few kpc, so the scales that 
are used to obtain the velocity dispersion are comparable. 
To confirm that the SDSS velocity dispersions are con¬ 
sistent with those obtained for the relatively small num¬ 
ber of individual galaxies with dynamically measured BH 
mass, we cross -correlated our pa r ent ga laxy sample with 
the catalog of lMcConnell fe Mai (1201311 . We found that 
three galaxies are listed in both samples. These values 
agree with each other within ~ 15%, thereby hinting 
that the central stellar velocity dispersions provided by 
SDSS can be used to trace the BH mass. Additionally, 
the tight scaling relation between the central stellar ve¬ 
locity dispersion and the BH mass is not only supported 
by empirical scaling relations, but both analytical and 
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hydrodynamical studies point out that this relation is 
a nat ural consequence of the self - regulated growth of 
BHs dSilk fe Reeslll998HKing|l2003t iSpringel et al.1 120051 

iHonkins et alJl2006HBooth fc Schavell2009ih 

Thus, both observational and theoretical considera¬ 
tions confirm that the central stellar velocity dispersion 
can robustly trace the average BH mass for the ellipti¬ 
cal galaxies in our study. To convert the central stel¬ 
lar velocity dispersion to BH mass, we used the best- 
fit Mrh — cr r relatio n obtained for elliptical galaxies 
(iGiiltekin et al.1120091 ). In Figure[5]we show the obtained 
scaling relation between the inferred BH and inferred to¬ 
tal gravitating mass within the virial radius (left panel) 
and within five effective radii (right panel). Addition¬ 
ally, in Figure [5] we also overplot measurements obtained 
for a sample of early -type galaxies with dyn amical BH 
mass measurements ((McConnelL&^M^JZOlSl) and accu¬ 
rately measured M (jKim fc Fabbianoll2013fl and M to t 
(Humphrey et al.1l2006l). 

To further investigate the correlation between BH and 
total halo masses, we compute the slope of the best-fit 
relation by implementing a least-squares (LS) linear re¬ 
gression to the data. Because simple LS regression in¬ 
cludes only the errors for the BH mass, we treat the halo 
mass for each data point as a normal distribution with 
width equal to the ler uncertainty shown in the insets of 
Figure [5] We randomly select a halo mass for each point 
from these distributions, and repeatedly refit the data 
using least-squares 5 x 10 4 times. The obtained best-fit 
relation between the inferred BH mass and the inferred 
total mass within the virial radius can be described as 

1 g + 0.6 

Mbh oc M tot 0,4 . The best-fit relation between the in¬ 
ferred BH mass and the inferred total mass within 5r e ff 

1 8 +0 ' 7 

can be expressed as Mbh cx M 5re “ 0 ' 6 . 

These results are in good agreement with previ¬ 
ous observat ional and theoretical studies. Indeed, 
iBandara et al.1 ( 2009 1 utilized galaxy-scale strong gravi¬ 
tational lenses to measure the total halo mass of galax¬ 
ies and observed Mr h oc M^f 5±031 . Moreov er, cosmo¬ 
logical simulations by iBooth fe Schavel (120101) predicted 
Mbh oc M^j) 5470 05 , and the authors suggested that the 
halo binding energy is the fundamental property that 
controls the BH mass. Overall, our results are consistent 
with this picture. Indeed, it is feasible that the central 
stellar velocity dispersion acts as an intermediary to di¬ 
rectly set the BH mass. However, as discussed in Section 
I5J1 the central stellar velocity dispersion is determined 
by the total halo mass. Thus, in elliptical galaxies the BH 
mass may be (indirectly) set by the ratio of the BH and 
the total halo masses. As opposed to these, the scaling 
relation observed between the BH and the stellar mass 
may b e the consequence of a velocity dispersion-to-halo 
mass (Iva n Uitert et alJl2 011h or the stellar-to -halo mass 
(|Mandelbaum et all 120061 : iMoster et all 1 20 1(1) relation. 

5.3. A Weaker Correlation between the BH and Stellar 

Masses 

To investigate the importance of the stellar mass in 
producing the observed X-ray luminosity, we repeat our 
stacking analysis by defining the galaxy bins in an in¬ 
verted manner. Specifically, we divide the galaxy sample 
into four equal-width central velocity dispersion bins in 



Fig. 6. — Results of the stacking analysis showing the stellar 
mass as a function of the X-ray luminosity for bins with fixed cen¬ 
tral velocity dispersion. The applied velocity dispersion ranges are 
marked with different symbols. The Mb u i ge — f-x relation reveals 
a weak correlation with a Pearson correlation coefficient of 0.73. 
Note that galaxies with similar central velocity dispersions have 
very similar X-ray luminosities despite their significantly different 
stellar masses. Along with the results of Figure^] this implies that 
the primary scaling relation is that between the BH mass and the 
total halo mass, and the stellar bulge plays a secondary role. 

the 110 — 310 km s _1 range. Similarly to our main analy¬ 
sis, these bins are further split based on their stellar mass, 
and the X-ray photons were stacked along with the cor¬ 
responding exposure maps. Following our earlier proce¬ 
dure, we extracted the source counts using a circular re¬ 
gion with four-pixel radius, and obtained the background 
level from a circular annulus with 8—14 pixel radii. For 
bins with X-ray luminosity of logLx/(erg s _1 ) > 40.1 
statistically significant (3.2er — 6.0cr) detections are ob¬ 
served, and for lower X-ray luminosities, moderately sig¬ 
nificant detections or nondetections are observed. The 
conversion of observed net counts to flux and the sub¬ 
traction of LMXBs were also performed as before (see 
Section 14.211 . 

Figure [6] shows the resulting Mb u i ge — Lx scaling re¬ 
lation, which reveals a significantly weaker correlation 
than that obtained between the central velocity disper¬ 
sion and the X-ray luminosity. Indeed, for a given stellar 
mass, the X-ray luminosity exhibits large variations: for 
example, at log M*/Mg ~ 10.9 we observe an order of 
magnitude difference in X-ray luminosity. We stress that 
this result refers to average galaxy populations because 
our stacking analysis, by definition, removes the distort¬ 
ing effects of individual outlier galaxies. The statistical 
analysis of the Mb u i ge — Lx relation produces a Pearson 
correlation coefficient of 0.73. For five degrees of freedom, 
the critical value to rule out Hq (i.e. no correlation be¬ 
tween these quantities) at only the 90% confidence level 
would require Hq > 0.81. Additionally, the best-fit rela¬ 
tion has a reduced % 2 = 2.2 for five degrees of freedom, 
which also hints at the large scatter. These results are 
in accord with earlier studies, which showed that for in- 
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dividual galaxies, at a given stellar mass, the luminosity 
of the hot X-ray gas may vary a s much as two orders of 
magnitude (iForman et al.l 119851: lO’Sullivan et al.1120011 : 

iMathews et al .1120031: iBoadan fc Gilfanovll2011h . Bevond 
the large scatter in the Mb u i g e — Tx relation, Figure [G] 
also reveals that galaxies with similar mean central stel¬ 
lar velocity dispersions have very similar X-ray luminosi¬ 
ties despite their markedly different stellar masses. This 
result is in excellent agreement with the conclusions pre¬ 
sented in Figure [H and illustrates that the stellar mass 
may play a more secondary role in determining the BH 
mass. 

In the present work we cannot probe whether a fun¬ 
damental plane exists that folds in the stellar mass of 
the elliptical galaxies. In this picture, the physical re¬ 
lation is dominated by the connection between the cen¬ 
tral potential and the dark matter, but also with the 
stellar mass play i ng a further, possibly secondary, role 
(lCapnellari et al.l 1201311 . However, given the shallow 
RASS observations, the connection between the inferred 
BH mass and total halo mass could only be studied in 
a relatively small number of bins. Because of this lim¬ 
itation, the further study of the fundamental plane is 
beyond the scope of the current paper. However, in the 
near future, the sensitive all-sky survey performed by 
eROSITA will allow us to split our galaxy sample into 
significantly smaller bins, and hence explore the param¬ 
eter space in greater detail. 

6. CONCLUSIONS 

In this work we utilized a sample of 3130 ellipti¬ 
cal galaxies observed by the SDSS and the RASS, and 
probed whether BHs exhibit a correlation with their dark 
matter halos. Our results can be summarized as follows: 

1. We grouped our galaxy sample based on stellar 
mass and central velocity dispersion, and found 
that the stacked X-ray images exhibit statisti¬ 
cally significant detections for bins with a c > 
190 km s” 1 and moderate detections or nondetec¬ 
tions for bins with cr c < 190 km s , irrespective of 
the median stellar mass of the galaxies in the bins. 
Based on the SDSS and stacked RASS data, we 
established that for massive elliptical galaxies the 
central stellar velocity dispersion exhibits a signifi¬ 
cantly tighter correlation with the X-ray luminosity 
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of the hot gas than with the stellar mass. 

2. We used the X-ray luminosity as a surrogate of the 
total halo mass, and concluded that it may be the 
central stellar velocity dispersion, and hence the 
central gravitational potential, that is most closely 
related to the dark matter halos of elliptical galax¬ 
ies. This hints that for elliptical galaxies the central 
stellar velocity dispersion may trace the large-scale 
dark matter halo more accurately than the stellar 
mass. 

3. By utilizing the central stellar velocity dispersion 
as a tracer of the BH mass, we established that the 
inferred BH and total halo masses tightly correlate 
with each other. The best-fit slopes between the in¬ 
ferred BH mass and the total mass within the virial 

i 6+°-6 

radius or 5 r e g can^ described as Mbh oc M tot 0,4 
or Mbh oc M 5re “ 0,6 , respectively. These results 
hint that in elliptical galaxies the central stellar 
velocity dispersion may act as an intermediary and 
directly set the BH mass. Because the central stel¬ 
lar velocity dispersion is determined by the dark 
matter halo mass, the BH mass may be (indirectly) 
set by dark matter halo mass. 
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